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Abstract 

We consider wireless networks that can be modeled by multiple access channels in which all the terminals 
are equipped with multiple antennas. The propagation model used to account for the effects of transmit and 
receive antenna correlations is the unitary-invariant-unitary model, which is one of the most general models 
available in the literature. In this context, we introduce and analyze two resource allocation games. In both 
games, the mobile stations selfishly choose their power allocation policies in order to maximize their individual 
uplink transmission rates; in particular they can ignore some specified centralized policies. In the first game 
considered, the base station implements successive interference cancellation (SIC) and each mobile station 
chooses his best space-time power allocation scheme; here, a coordination mechanism is used to indicate to 
the users the order in which the receiver applies SIC. In the second framework, the base station is assumed to 
implement single-user decoding. For these two games a thorough analysis of the Nash equilibrium is provided: 
the existence and uniqueness issues are addressed; the corresponding power allocation policies are determined 
by exploiting random matrix theory; the sum-rate efficiency of the equilibrium is studied analytically in the low 
and high signal-to-noise ratio regimes and by simulations in more typical scenarios. Simulations show that, in 
particular, the sum-rate efficiency is high for the type of systems investigated and the performance loss due to 
the use of the proposed suboptimum coordination mechanism is very small. 
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I. Introduction 

In this paper, we consider the uplink of a decentralized network of several mobile stations (MS) and one base 
station (BS). This type of network is commonly referred to as the decentralized multiple access channel (MAC). 
The network is said to be decentralized in the sense that each user can freely choose his power allocation (PA) 
policy in order to selfishly maximize a certain individual performance criterion, which is called utility or payoff. 
This means that, even if the the BS broadcasts some specified policies, every user is free to ignore the policy 
intended for him if the latter does not maximize his performance criterion. 

To the best of the authors' knowledge, the problem of decentralized PA in wireless networks has been 
properly formalized for the first time in HI, (H. Interestingly, this problem can be formulated quite naturally 
as a non-cooperative game with different performance criteria (utilities) such as the carrier-to-interference ratio 
0, aggregate throughput [|4l or energy efficiency |0, O. In this paper, we assume that the users want to 
maximize information-theoretic utilities and more precisely their Shannon transmission rates. Indeed, the point 
of view adopted here is close to the one proposed by the authors of Q for DSL (digital subscriber lines) 
systems, which are modeled as a parallel interference channel; (H for the single input single output (SISO) 
and single input multiple output (SIMO) fast fading MACs with global CSIR and global CSIT (Channel State 
Information at the Receiver/Transmitters); Q for MIMO (Multiple Input Multiple Output) MACs with global 
CSIR, channel distribution information at the transmitters (global CDIT) and single-user decoding (SUD) at 
the receivers; iflOl . ifTTI for Gaussian MIMO interference channels with global CSIR and local CSIT and, by 
definition of the conventional interference channel lfT2l . SUD at the receivers. Note that reference |fT3l where 
the authors considered Gaussian MIMO MACs with neither CSIT nor CDIT differs from our approach and that 
of 0, JH, (9l, (TOl, Ell because in HH the MIMO MAC is seen as a two-player zero-sum game where the 
first player is the group of transmitters and the second player is the set of MIMO sub-channels. The closest 
works to the work presented here are [9ij and |[T4l . Although this paper is in part based on these works, it still 
provides significant contributions w.r.t. to them, as explained below. 

In (91, the authors consider MIMO multiple access channels and assume SUD at the BS; the authors formulate 
the PA problem into a team game in which each user chooses his PA to maximize the network sum-rate. In 
lfl4l . the same type of decentralized networks is considered but SIC is assumed at the BS. As each user needs 
to know his decoding rank in order to adapt his PA policy to maximize his individual transmission rate, a 
coordination mechanism has to be introduced: the coordination signal precisely indicates to all the users the 
decoding order used by the receiver. The present paper differs from these two contributions on at least four 
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important technical points: (i) when SUD is assumed, the PA game is not formulated as a team game but as 
a non-cooperative one; (ii) we exploit several proof techniques that are different from @; (iii) while and 
Ifl4l assume a Kronecker propagation model with common receive correlation we assume here a more general 
model, the unitary-invariant-unitary (UIU) propagation model introduced by 11211 . for which the users can have 
different receive antenna correlation profiles. This is useful in practice since, for instance, it allows one to 
study propagation scenarios where some users can be in line of sight with the BS (the receive antenna are 
strongly correlated) whereas other users can be surrounded by many obstacles, which can strongly decorrelate 
the receive antennas for these users; (iv) while the authors of Ifl4l restricted their attention to either a purely 
spatial PA problem or a purely temporal PA problem, we tackle here the general space-time PA problem. 

In this context, our main objective is to study the equilibrium of two power allocation games associated with 
the two types of decoding schemes aforementioned (namely SIC and SUD). The motivation for this is that the 
existence of an equilibrium allows network designers to predict, with a certain degree of stability, the effective 
operating state(s) of the network. Clearly, in our context, uniqueness is a desirable feature of the equilibrium. 
As it will be seen, it is possible to prove the existence in both games under investigation. Uniqueness is 
proven in the case of SUD while it is conjectured for the case of SIC. In order to establish the corresponding 
results, the paper is structured as follows. After presenting the general system model in Sec. J]] we analyze 
in detail the space-time PA game when SIC and a corresponding coordination mechanism are assumed (Sec. 
ITTll ). For this game, the existence and uniqueness of the NE are proven and the equilibrium is determined by 
exploiting random matrix theory when the numbers of antennas are sufficiently large. Its sum-rate efficiency 
is also analyzed. In Sec. [TV] we analyze the case of SUD since this decoding scheme, although suboptimal in 
terms of performance (even in the case of a network with single-antenna terminals), has some features that can 
be found desirable in some contexts: the receiver complexity is low, there is no need for a coordination signal, 
there is no propagation error since the data flows are decoded in parallel and not successively and also it is 
intrinsically fair. To analyze the case of the SUD-based PA game, we will follow the same steps as in Sec. [TTT] 
and we will see that, the equilibrium analysis can be deduced, to a large extent, from the SIC case. Numerical 
results are provided in Sec. [V] to illustrate our theoretical analysis and to better assess the sum-rate efficiency 
of the considered games. Sec. [Vl] corresponds to the conclusion. 

II. System Model 

We assume a MAC with arbitrary number of users, K > 2. Regarding the original definition of the MAC by 
lfl31 and lfl6l . the system under consideration has two common features: all transmitters send at once and at 
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different rates over the entire bandwidth, and the transmitters are using good codes in the sense of the Shannon 
rate. Our system differs from |[I31l[l6l in the sense that multiple antennas are considered at the terminal nodes, 
channels vary over time and the BS does not dictate the PA policies to the MSs. Also, we assume the existence 
of coordination signal which is perfectly known to all the terminals. If the coordination signal is generated by the 
BS itself, this induces a certain cost in terms of downlink signaling but the distribution of the coordination signal 
can then be optimized. On the other hand, if the coordination signal comes from an external source, e.g., an FM 
transmitter, the MSs can acquire their coordination signal for free in terms of downlink signaling. However this 
generally involves a certain sub-optimality in terms of uplink rate. In both cases, the coordination signal will be 
represented by a random variable denoted by S G S. Since we study the K— user MAC, S = {0, 1, Kl} is 
a K\ + 1-element alphabet. When the realization is in {1, K\} , the BS applies SIC with a certain decoding 
order (game 1). When S = the BS always applies SUD (game 2), where all users are decoded simultaneously 
(no interference cancellation). In a real wireless system the frequency at which the realizations would be drawn 
would be roughly proportional to the reciprocal of the channel coherence time (i.e., l/T co h). Note that the 
proposed coordination mechanism is suboptimal because it does not depend on the realizations of the channel 
matrices. We will see that the corresponding performance loss is in fact very small. 

We will further consider that each mobile station is equipped with n t antennas whereas the base station has 
n r antennas (thus we assume the same number of transmitting antennas for all the users). In our analysis, the 
flat fading channel matrices of the different links vary from symbol vector (or space-time codeword) to symbol 
vector. We assume that the receiver knows all the channel matrices (CSIR) whereas each transmitter has only 
access to the statistics of the different channels (CDIT). The equivalent baseband signal received by the base 
station can be written as: 

K 

y«(r) = ^H k (r)xi s \T) + Z^(r), (1) 

k=l 

where (r) is the n t -dimensional column vector of symbols transmitted by user k at time r for the realization 
s G S of the coordination signal, H^(r) € C™ rXnt is the channel matrix (stationary and ergodic process) of 
user k and Z^ s \t) is a n r -dimensional complex white Gaussian noise distributed as N(0, a 2 I nr ). For the sake 
of clarity we will omit the time index r from our notations. 

In order to take into account the antenna correlation effects at the transmitters and receiver, we will assume 
the different channel matrices to be structured according to the unitary-independent-unitary model introduced 
in (HI: 

VA; e {!,..., K}, H k = V k H k W k , (2) 
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where and Wfc are deterministic unitary matrices that allow one to take into consideration the correlation 
effects at the receiver and transmitter. Also is an n r x nt matrix whose entries are zero-mean independent 
complex Gaussian random variables with an arbitrary profile of variances, such that E\H k (i, j)| 2 = 2*M The 



Kronecker propagation model for which the channel transfer matrices factorizes as 



Rl /2 ®k:T 



1/2 



is a 



,2 _ 4 R) (*)4 T) Q-) 



special case of the UIU model where the profile of variances is separable i.e., E\H k (i,j)\' z = " k ^ u ' , with 
for each k: @k is a random matrix with zero-mean i.i.d. entries, is the transmit antenna correlation matrix, 
Rfc is the receive antenna correlation matrix, {d k (j)}je{i,...,n t } an ^ {% W}ie{i,-,n r } are tne i r associated 
eigenvalues. In this paper we will consider that V^, = V for all users. The reason for assuming this will be made 
clearer a little further. In spite of this simplification, we will still be able to deal with some useful scenarios 
where the users see different propagation conditions in terms of receive antenna correlation. 



III. Successive Interference Cancellation 
When SIC is assumed at the BS, the strategy of user k € {1, 2, K}, consists in choosing the best vector 



x (s) x (s),H 
— k — k 



, for s € S, in the sense of 



of precoding matrices = ^QJ,, ,Q k , Q k ) where Q k — E 
his utility function. For clarity sake, we will introduce another notation which will be used in the remaining 
of this section to replace the realization s of the coordination signal. We denote by Vk the set of all possible 
permutations of K elements, such that ir € Vk denotes a certain decoding order for the K users and ir(k) 
denotes the rank of user k G K, and it~ 1 € Vk denotes the inverse permutation (i.e. 7r~ 1 (7r(/c)) = k) such that 
7r~ 1 (r) denotes the index of the user that is decoded with rank r € /C. We denote by p n S [0, 1] the probability 
that the receiver implements the decoding order ir £ Vk, which means that p n = 1. At last note that there 

ttEVk 

is a one-to-one mapping between the set of realizations of the coordination signal S and the set of permutations 
Vk, i-e. £ : S — > Vk such that £(•) is a bijective function. This is the reason why the index s can be replaced 
with the index n without introducing any ambiguity or loss of generality. The vector of precoding matrices can 
be denoted by Q = ( ) and the utility function can be written as: 



where 



..sic 
u k 



(Qk,Q-k)= £ rf } (Q^Q-t 



Elog 2 



I + pH fc Q«Hf + p £ H.Q^H 



Mi 



Elog 2 



I + P 



E 



(3) 



(4) 
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with p = \ and K, k = £ )C\ir(£) > 7r(k)} represents, for a given decoding order n, the subset of users 
that will be decoded after user k. Also, we use the standard notation —k , which stands for the other players 
than k. An important point to mention here is the power constraint under which the utilities are maximized. 
Indeed for user k £ {1, ...,K}, the strategy set is defined as follows: 

Af c = lq k =(Q ( ! : ) )^ v | VneV K ,Q { k ] hO, ^ ^Tr(Q«) < n t P fc | . (5) 

tt£V k ) 

In order to tackle the existence and uniqueness issues for Nash equilibria in the general space-time PA game, 
we exploit and extend the results from Rosen ifTTl , which we will briefly state here below in order to make 
this paper sufficiently self-contained. 

Theorem I: IfTTl Let Q = (fC, {Ak}keic, {^fc}fce/c) ^ e a game where K = {1,...,K} is the set of players, 
Ax, ...,Ak the corresponding sets of strategies and u\, ...,Uk the utilities of the different players. If the following 
three conditions are satisfied: ( i) each u k is continuous in the all the strategies aj £ Aj , Vj £ /C; ( ii) each u k 
is concave in a k € A k ; (Hi) A\, ...,Ak are compact and convex sets; then Q has at least one NE. 

Theorem 2: M17II Consider the K -player concave game of Theorem [7J If the following (diagonally strict 

concavity) condition is met: for all k £ K, and for all (af k ,a' k ) £ A\ such that there exists at least one index 

K 

j £ K for which (J- ^ d-, — af k ) T [V 'a k Uk{of k , d~k) ~ ^ a k u k{s^i fi-fe)] > 0> men me uniqueness of the 

k=l 

NE is insured. 

In the space-time power allocation game under investigation, the obtained results are stated in the following 
theorem. 

Theorem 3: [Existence of an NE] The joint space-time power allocation game described by: the set of players 
k £ K,; the sets of actions A k and the utility functions u k lc (Q k , Q-fc) given in (0), has a Nash equilibrium. 

Proof: It is quite easy to prove that the strategy sets A k are convex and compact sets and that the utility 
functions u k lc (Q k , Q-k) are concave w.r.t. Q k and continuous w.r.t. to (Q^, Q-k) and by Theorem[TJat least 
one Nash equilibrium exists. For more details, the reader is referred to Appendix |A] ■ 

Theorem 4: [ Sufficient condition for uniqueness ] If the following condition is met 

E E Tr " ( v q^ ic (<^ QU) - v Q (^ IC (Q£, Q-k)) } > o (6) 

tt£Vk k=l 

for all Q' k = (Q k nY ) , Q'l = (qP") £ Af c such that (Q[, . . . , Q' K ) ^ (Q'(, Q^), then the 
Nash equilibrium in the power allocation game of Theorem \3\ is unique. 

This theorem corresponds to the matrix generalization of the diagonally strict concavity (DSC) condition of 
[17] and is proven in Appendix B. To know whether this condition is verified or not in the MIMO MAC one 
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needs to re-write it in a more exploitable manner. It can be checked that C expresses as C = P-kT-k where 
for each it € Vk, T-k is given by: 

% = £Tr{(qr-Qr) ^^(qw'.qw 1 )-v Qr ,i?r ) (Qr,QLT)]} 
fe=i 

= e j> { P H^ w (Qw;; r) - Q^; (r) )Hf_ 1(r) 

(i+pH^ 1(r) Qw; (r) Hf_ 1(r) + P f; vw^Hf-J 

V s=r+l J 

I + P H »r-Hr)Q ff - ) i(p) H ff-i(r) + P XI H t 1 WQi-i^) 11 ^ 1 (s) ) 

V s=r+l / 



(7) 



A* 



eJ^thaW" -aW 



r=l 



A' 



A 



where A 



to' 



= E[^(H)] 

pH 7r -i (r) Q^ ) 1(r) H^_ 1(r , ) , A^" ; = pH^-i^Q^^H^L^ and the users have been ordered 
using their decoding rank rather than their index. Notice that since the expectation operator is linear we can 
switch between the trace and expectation. 



Let us denote by H = [H 1 ,...,Ha], Q' = (Q' k ) kelc , , Q" = (Q't) k&c , Q W ' = (QF ), r > ' Q 



fce/c 



, A 



W 



n£V K ,kelC 



, A" 



(tt)" 



fce/c 



tts-Pk, fce/c 



In order to prove that the DSC condition holds we have to prove that for all Q' ^ Q" we have C > 0. 
Let us give a very useful result. 

Lemma 1: For any positive definite matrices Ai, Bi, ant/ any positive semi-definite matrices Aj, B^, 
i € {2, . . . , K}, we /zave f/iaf 

K 



^ Tr < (Ai — Bj 



i=l 



5> 



> 



(8) 



where the equality holds if and only if Aj = B^ for all j € {1, ... , K} 

The proof can be found in ll27ll . for K = 2, and in |[28l for arbitrary K > 2. Using this result, we can prove 
that for any channel realization, any Q', Q" and any 7r € Vk- 



F n (U) = Tr ( A«" - AW' 



> 



(9) 
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implying that 7^ > and that C > 0. Let us consider now two arbitrary covariance matrices such that Q' ^ Q". 
This means that there is at least one decoding order $ € Vk such that QW ^ QW". We will prove that 
T$ > which will imply the desired result C > 0. 

Remark: Assuming that rank(H^Hfc) = n t , for all k £ /C, and n< < n r + 1, then Q' ^ Q" implies that 
A' 7^ A". This means that for any channel realization we have -F$(H) > which implies directly Tq > and 
C > 0. 

For the general proof, let us define the following sets: 

A H (Q W ,Q W ") = {H€V H \\/k€K:K k (Q k W - Q k W)H^ = o} 



(10) 



We know that: 



% = E[F„(H)] 

= / Fa(H)L(H)dH 

= / F^(H)L(H)dH+ / i^(H)L(H)dH 

■^^ff(Q""',Q (< "") ■AAff(Q w ',QW") 

F^(H)L(H)dH 



(11) 



U H (QW',QW") 

where L(H) > stands for the p.d.f. of H € T>h = C nrXKnt . The second equality follows since T>h = 
A H (Q W ', Q w ')UAh(Q W , Q W ")- The third equality follows since for all H € A H (Q W , Q W ") we have 
that Ftf(H) = from Lemma [8] We know that for all H e ^ H (Q W ', Q W ") we have that ify(H) > 0. It 
suffices to prove that A H (Q W ',Q W ") is a subset of non-zero Lebesgue measure to imply that T$ > and 
thus that C > 0. It turns out that we can prove the existence of a compact set C Ah (Q^ > Q^ ") for 
arbitrary QW ^ QW". Thus, we have the desired result C > 0. 

Determination of the Nash equilibrium. In order to find the optimal covariance matrices, we proceed in the 
same way as described in First we will focus on the optimal eigenvectors and then we will determine the 
optimal eigenvalues by approximating the utility functions under the large system assumption. 

Theorem 5: [ Optimal eigenvectors ] For all k € K,, Q k € Af lG there is no loss of optimality by imposing 
the structure Q k = (Qk^)iv£VK> Qfc^ = W k P k ^ W k H , in the sense that: 

max ^ IC (Q fc , Q_ fc ) = max ^ IC (Q fc , Q_ k ), 



where Sf c = {q fe = (Q k (w) )«&> h G Af c \Q k n) = W fc P^ } Wf }, s € S, model from (TJ) and p[ s) 
Diag(pW(l),... ) P fe W K)). 
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The detailed proof of this result is given in Appendix This result, although easy to obtain, it is instrumental 
in our context for two reasons. First, the search of the optimum precoding matrices boils down to the search of 
the eigenvalues of these matrices. Second, as the optimum eigenvectors are known, available results in random 
matrix theory can be exploited to find an accurate approximation of these eigenvalues. Indeed, the eigenvalues 
are not easy to find in the finite setting. They might be found using numerical techniques based on extensive 
search. Here, our approach consists in approximating the utilities in order to obtain expressions which are not 
only easier to interpret but also easier to be optimized w.r.t. the eigenvalues of the precoding matrices. The key 
idea is to approximate the different transmission rates by their large-system equivalent in the regime of large 
number of antennas. The corresponding approximates can be found to be accurate even for relatively small 
number of antennas (see e.g., lfT8llfT9l for more details). 

Since we have assumed V*. = V, we can exploit the results in EOlETl for single-user MIMO channels, 
assuming the asymptotic regime in terms of the number of antennas: n r — > oo, n t — > oo, ^ — > (3. The 
corresponding approximated utility for user k is: 




(12) 



where 





(13) 



n 
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where = |/c£° | and the parameters 7™ (j) and 5^> (j) Vj 6 {1, . . . , n t }, k G K, tt G Vk are the solutions 
of: 

yje{l,...,n t },£eJC { ^U{k}: 

1 ^ cr e (i,j) 



?(*")/ 



7f(i) 



'"^+l)n t 



E 



(ivW + l)pP^(j) 



+ i)„ t ^ ^ ffr (i,m)5W(m) 
re/c^'ulfc}^ 1 



(14) 



i+(N^+i) P pr(jhr(j) 



and <j>iij), Vj G {1, . . . ,n t } and 7r G Pjf are the unique solutions of the following system: 



VjGil,...,^},^^ = 



(j) 



M m =l 



The corresponding water-filling solution is: 



(15) 



p W,N E(i) 



1 



1 



ln2n r A fc pj^ ti) 

where A& > is the Lagrangian multiplier tuned in order to meet the power constraint: 

1 1 



n t P k . 



(16) 



ln2n r A fc N^py^ti) 

Note that to solve the system of equations given above, we can use the same iterative power allocation algorithm 
as the one described in |9). 

At this point, an important point has to be mentioned. The existence and uniqueness issues have be analyzed 
in the finite setting (exact game) whereas the determination of the NE is performed in the asymptotic regime 
(approximated game). It turns out that large system approximates of ergodic transmission rates have the same 
properties as their exact counterparts, as shown recently by |f23l , which therefore ensures the existence and 
uniqueness of the NE in the approximated game. 

Nash Equilibrium efficiency. In order to measure the efficiency of the decentralized network w.r.t. its cen- 
tralized counterpart we introduce the following quantity: 



dNE 

SRE = ±» < 1, 



(17) 
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where SRE stands for sum-rate efficiency; the quantity represents the sum-rate of the decentralized network 
at the Nash equilibrium, which is achieved for certain choices of coding and decoding strategies; the quantity 
Csum corresponds to the sum-capacity of the centralized network, which is reached only if the optimum coding 
and decoding schemes are known. Note that this is the case for the MAC but not for other channels like the 
interference channel. Obviously, the efficiency measure we introduce here is strongly connected to the price 
of anarchy ll24l (POA). The difference between SRE and POA is subtle. In our context, information theory 
provides us with fundamental physical limits on the social welfare (network sum-capacity) while in general no 
such upper bound is available. In our case, the sum-capacity is given by: 



C sum = max E log 
(n 1 ,...,n K )eA^ 



K 



i + P £)H fc n fc H£ 



fe=i 



(18) 



with 

A {c) = {(ni,...,nji-)|v* g ic,n k h o,n k = n%,Tr(n k ) < n t p k } . (19) 

In general, it is not easy to find a closed-form expression of the SRE. This is why we will respectively analyze 
the SRE in the regimes of high and low signal-to-noise ratio (SNR), and for intermediate regimes simulations 
will complete our analysis. It turns out that the SRE tends to 1 in the two mentioned extreme regimes, which 
is the purpose of what follows. 

In the high SNR regime, where p — s- oo, we observe from (fT4l ) that 8^\j) — > {7r ) . Under this condition, it 
is easy to check that by setting the derivatives of C k w.r.t. Pj, s \j) to zero, we obtain that the power allocation 
policy at the NE is the uniform power allocation p(' r )' NE = p k j t regardless the realization of the coordination 
signal S. Furthermore, in the high SNR regime, the sum-capacity is achieved by the uniform power allocation. 
Thus, we obtain that the gap between the NE achievable sum-rate and the sum-capacity is optimal, SRE = 1 
for any distribution of S. 

In the low SNR regime, where p — s> 0, from (fl4l) we obtain that 5^ (j) — > and that 7^ (j) = ^ — ae(i, j). 

i=l 

By approximating ln(l + x) x when x « 1, the power allocations policies at the NE are the solutions of 



the following linear programs: 



{P k te)h<i<nt j=1 UeP K i=i 



(20) 



s -t- E E Pk\j)<P k nt 
j=i Trev K 
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given by: 



£ P^ w,NE (i) 



?i t P fc ifj = arg max V" a k (i, m) 

l<m<n t *■ — ' 
i=l 



(21) 







otherwise 



The optimal power allocation that achieves the sum-capacity is equal to the equilibrium power allocation, 
P* k = YlireV Pn*k ^ NE (i) Thus, the achievable sum-rate at the NE is equal to the centralized upper bound 
and thus SRE = 1 for any distribution of S. In conclusion, when either the low or high SNR regime is 
assumed, the sum-capacity of the fast fading MAC is achieved at the NE although a sub-optimum coordination 
mechanism is assumed and also regardless of the distribution of the coordination channel. 



IV. Single User Decoding 

In this section the coordination signal is deterministic (namely Pr[5 = s] = 5(s), 5 being the Kronecker 
symbol) and therefore the amount of downlink signalling the BS needs in order to indicate to the MSs that it is 
using SUD can be made arbitrary small (by letting the frequency at which the realizations of the coordination 
signal are drawn tend to zero). In this framework, each user has to optimize only one precoding matrix. Indeed, 



the strategy of user k G /C, consists in choosing the best precoding matrix Q 



(o) 



E 



X (0) X (0)H 



, in the sense 



of his utility function obtained with SUD: 



^ UD (Qi 0) ,Q^) = Ei og 



I + pH fc Qi°>H? +p£HtQ? ) B? 
. The strategy set of user k becomes 



Elog 



A f» = | Q (0) ^ O Q (0) = Q (0)^ )Tr(Q (0)) < nt p fc } . 



(22) 



(23) 



It turns out that the equilibrium analysis in the game with SUD can be, to a large extent, deduced from the 
game with SIC. For this reason, we will not detail the corresponding proofs. The existence and uniqueness 
issues are given in the following theorem. 

Theorem 6: [Existence and uniqueness of an NE] The space power allocation game described by: the set 
of players k € 1C; the sets of actions -4| UD and the payoff functions ^fc UD (Q[°\ Q-l) given in ri22D , has a 
unique Nash equilibrium. 

To prove the existence of a Nash equilibrium we also exploit Theorem [TJ and the four necessary conditions 
on the utility functions and strategy sets can be verified using the same tools as described in Appendix |A] 

Uniqueness of the Nash equilibrium. Here we can specialize Theorem 4, which is the matrix extension of 
Theorem 2. When the strategies sets are not sets of pairs of matrices but only sets of matrices, the diagonally strict 
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concavity condition in (6) can be written as follows. For all q[ 0) ', Q^ 0) " G Al VB such that (qJ 0) ', . . . , Q^') ^ 



(QS 0) ",...,Qf': 



fc=i 

Now we can evaluate C and obtain that: 



fc=i 



C = E^Tr I P H fc (Q( 0) '-Qf")H 



K 



I + pJ> Q f" H f 




(25) 



A 



= ETr{(B'-B")[(B")- 1 -(B')- 1 ]}, 
= E[F (H)] 

if 

which is positive for any B = I + ^H^Q (0) 'Hf , B = I + ^H f Q (0) "Hf from © for K = 2. We need 

i=\ t=\ 
to prove that for any (Cjf ) ', . . . , Q^') ^ (qJ° } ", . . . , Q^") we have C > 0. 

Remark: Assuming that rank(H^H) = Kn t and ifjij <n r + K, then Q' ^ Q" implies that B' ^ B". 

This means that for any channel realization we have -Fo(H) > which implies directly that C > 0. 

For the general proof, we define the following sets: 

MQ (0) ',Q (0) ") = {He^|Ef =1 H fc (Q fc (°)'-Q fc (°)")H^ = o} 

MQ (0) ',Q (0) ") = {HG2? H |EjLiH fc (Q fc (°)'-Q fc (°)")H^ #0} 

We know that: 



(26) 



(27) 



C = E[F (H)] 

= / F (H)L(H)dH 

= I F (H)L(H)dH + / F (H)L(H)dH 

JB H (Q m ' ,Q m ") JB a W°y ,QC°)") 

= / F (H)L(H)dH 

^Bff(Q <0) ',Q (0) ") 

The second equality follows since V H = B H (Q {0) ' , Q {0) ") U S//(Q {0) ', Q (0) ")- The third equality follows 
because F (H) = for all H G Bff(Q (0) ',Q (0) ") from Lemma 1. We also know that F (H) > for all 
H G B H (Q {oy , Q (0) "). It suffices to prove that B H (Q {oy , Q (0) ") is a subset of non-zero Lebesgue measure to 
imply that C > 0. Here as well, the existence of the compact set can be proved (similarly to the proof for the 
SIC decoding technique). 
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Determination of the Nash equilibrium. As for the optimal eigenvectors of the covariance matrices, we follow 
the same lines as in Appendix In this case also there is no loss of optimality by choosing the covariance 
matrices = W^Pj^Wj^, where is the same unitary matrix as in Q and is the diagonal matrix 
containing the eigenvalues of q! . 

Here also we further exploit the asymptotic results for the MIMO channel given in EOl l2ll . The approximated 
utility for user k is: 



,~;SUD/p(0) p (0). 



K m 



- E E ^(i + K P pf\j) lk (j))+ 

r k=l j=l 

j n r / K n t \ 

-E lo g 2 h+^DE^w - 

t=i \ fc=i i=i / 

^ K nt 

—EE ik{j)h{j) lo §2 e ~ 

r fc=i j=i 
nt 

-EE ^(i + (* - VpH%)M))- 

T t+k 3=1 

1 n r I 1 n t \ 

-X> 2 i+^EE^^oi J + 



71 



i;EE^W*W log 2 e 
where the parameters 7fc(j) and 5fc(j) Vj € {1, ... , n t }, k E {1, 2} are solution of: 



Vje{l,...,nt},fce£: 



TfcC?) 



1 



£ 

i=i 



o-k(i,j) 



K rn 



1 + ^;E E Mh m ) s i( m ) 



=1 m=l 



l + ^pP, (0) (i) 7fc (j-) 



and (f>i(j), ipi(j), Vj € {1, . . . , nt} are the unique solutions of the following system: 



Vj G {l,...,ntye/C\W : 



i 



E 



o-e(i,j) 



(K - l)n t * 

1 + I^TkE E o-r(*,m)^ r (m) 

r^fc m=l 

l + (iC-l)pR (0) (j)^(i)' 



(28) 



(29) 



(30) 
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The corresponding water-filling solution is: 




ln2n r A fc Kpj k (j) 



1 1 



(31) 



+ 



In what the efficiency of the NE point is concerned, we already know that the SUD decoding technique 
is sub-optimal in the centralized case (SUD does allow the network to operate at an arbitrary point of the 
centralized MAC capacity region) and it is impossible to reach the sum-capacity C sum even if the high and low 
SNR regime are assumed. 



In what follows, we assume the regime of large numbers of antennas. From O, |[20l , (TO, we know that 
the approximates of the ergodic achievable rates in the asymptotic regime are accurate even for relatively small 

1/2 1/2 

number of antennas. For the channel matrices, we assume the Kronecker model = R fc & k T k mentioned 
in Sec.ini where the receive and transmit correlation matrices R^, follow an exponential profile characterized 



Tfc(i,i) = tl . By assuming that the receive antenna is a uniform linear array (ULA) and knowing that, when 
the dimensions of Toeplitz matrices increase they can be approximated by circular matrices we obtain that all 
the receive correlation matrices can be diagonalized in the same vector basis (i.e., the Fourier basis). Thus 
the considered model is included in the UIU model that we studied where = V. 

Fair SIC decoding versus SUD decoding. First we compare the results of the general space-time PA game 
considered in Sec. [Till where SIC decoding is used at the receiver, and the game described in Sec. UVj where 
SUD decoding is used. Fig. Q] depicts the achievable sum-rate at the equilibrium as a function of the transmit 
power Pi = Pi = P, for the scenario n r = n t = 10, r = [0.5, 0.2], t = [0.5, 0.2], p = 3dB. In order to have a 
fair comparison we assume that p = | (on average each user is decoded second half of the time when SIC is 
assumed). We observe that, even in this scenario, which was thought to be a bad one in terms of sub-optimality, 
the sum-rate obtained with the first game is very close to the sum-capacity upper bound. Also, the sum-rate 
reached when the BS uses SUD is clearly much lower than the sum-rate obtained by using SIC. 

SIC decoding, comparison between the joint space-time PA and the special cases of spatial PA and temporal 
PA. Now we want to compare the results of the general space-time PA with the two particular cases that were 



V. SIMULATION RESULTS 



by the correlation coefficients (see e.g., |[25l . |[26l ) r = [ri,r2] and t = [<i,<2] such that H k (i,j 
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studied in lfl4l : the spatial PA, where the users are forced to allocate their power uniformly over time (regardless 
of their decoding rank) but are free to allocate their power over the transmit antennas; the temporal PA, where 
the users are forced to allocate their power uniformly over their antennas but they can adjust their power as 
a function of the decoding rank at the receiver. Fig. [2] represents the sum-rate efficiency as a function of the 
coordination signal distribution parameter p € [0, 1] when n r = m = 10, r = [0.3,0], t = [0.5,0.2], p = AdB, 
Pi = 5, P2 = 50. We observe that the three types of power allocation policies perform very close to the upper 
bound. What is most interesting is the fact that the performance of the network at the equilibrium is better by 
using a purely spatial PA instead of the most general space-time PA. This has been confirmed by many other 
simulations and illustrates a Braess paradox: although the sets of strategies for the space-time case include 
those of the purely spatial case, the performance obtained at the NE are not better in the space-time case. 

SIC decoding, spatial PA, achievable rate region. In Fig. [3j we observe that the rate region achieved at the 
NE of the space PA as a function of the distribution of the coordination signal p for the scenario n r = n t = 10, 
r = [0.4,0.2], t = [0.6,0.3], p = MB, P 1 = 5, P 2 = 50. It is quite remarkable that in large MIMO MACs, the 
capacity region comprises a full cooperation segment just like the SISO MACs. The coordination signal precisely 
allows one to move along the corresponding line. This shows the relevance of large systems in decentralized 
networks since they allow to determine the capacity region of certain systems whereas it is unknown in the 
finite setting. Furthermore, they induce an averaging effect, which makes the users' behavior predictable. 

VI. CONCLUSIONS 

Interestingly, the existence and uniqueness of the Nash equilibrium can be proven in multiple access channels 
with multi-antenna terminals for a general propagation channel model (namely the unitary-invariant-unitary 
model) and the most general case of space-time power allocation schemes. In particular, the uniqueness proof 
requires a matrix generalization of the second theorem of Rosen flTTl and proving a trace inequality ||28l . For 
all the types of power allocation policies (purely temporal PA, purely spatial PA, space-time PA), the sum-rate 
efficiency of the decentralized network is close to one when SIC is assumed and the network is coordinated by 
the proposed suboptimum coordination mechanism. Quite surprisingly, the space-time power allocation performs 
a little worse than its purely spatial counterpart, which puts in evidence a Braess paradox in the types of wireless 
networks under consideration. One of the interesting extensions of this work would be to analyze the impact of 
a non-perfect SIC on the PA problem. Indeed, the effect of propagation errors could then be assessed (which 
does not exist with SUD). 
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Appendix A 



A. Concavity of the utility functions u k lc 



Let us focus on user k € /C. We want to prove that ^ IC (Qfc, Q-fc) is concave w.r.t. € ^1^ IC . We observe 
that the term (q! , Q_2) in © depends only on and and not on the covariance matrices Q k T \ 



Q-l f° r an y otner possible decoding rule r € Vk \ {^}- Thus, in order to prove that ii| IC (Qfc, Qk) is sttictly 



concave w.r.t. to = (Q k ) 7rG -p K , it suffices to prove that R k (Q k ,Q_D is concave w.r.t. Q k for all 

To this end, we study the concavity of the function /(A) = ^(AQ^' + (1 - A)Q^ ) ") over the interval 
[0, 1] for any pair of matrices (Q^ , )■ The second derivative of / is equal to: 



-ETr 



p 2 Hf I I + pH fc QW"Hf + pAHfcAQ^Hf + p £ H.Q^Hf I H fe AQ^ 



(f)lTil 



xHf I I + pH fc Qf "Hf + pAH fc AQf Hf + p £ H,Q«Hf ) H fe AQ« 
= -ETrfAAQ^AAQ^] 

with A = p 2 Hf [ I + pH A .Q^ ) "Hf + pAHjAQ^Hf + p ^ H^Q^Hf j H fc , which can be proven 



to be a Hermitian positive definite matrix, AQ fc 



— Q^. also a Hermitian matrix, and p = 



&2f (\) 



-ETr[A 1 / 2 AQ^A 1 / 2 A 1 / 2 AQ^A 1 / 2 ] 
-ETr[BB^] < 



with B = A 1 / 2 AQWA 1 /2. 



B. Continuity of the utility functions u k 

Considering the Leibniz formula, the determinant of a matrix can be expressed as a weighted sum of products 
of its entries. Knowing that the product and the sum of continuous functions are continuous, we conclude that 
the determinant function is continuous. Also, it is well known that the logarithmic function is a continuous 
function. Thus, for any it € Vk, the function i?i (qL , Q^) is nothing else but the composition of two 
continuous functions which is also continuous w.r.t. (Q k n \Q^ k ). This suffices to prove that u k lc (Qk, Q-k) 
is continuous w.r.t. (Qfc,Q_fc). 
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C. Convexity of the strategy sets *4| IC 

In order to prove that the set w4| IC is convex, we need to verify that, for any two matrices (Q fc ,Q^) G 
Af c x Af c , we have: 

aQ\ + (1 - a)Q'i G Af c , 

for all a > 0. 

For any Q fe ,Q^ G A^ lC \ the matrices Qj^ are Hermitian which implies that aQ^ + (1 — a>)Q^ are 
also Hermitian matrices, for all tt G Vk- 

Furthermore, for any Ql, G -4| IC , we have that , q[ are non-negative matrices which implies 
that aQ^ + (1 — a)Q^ are also non-negative matrices, for all 7r G Vk- 

Finally, knowing that the trace is a linear application we have that: 

J>.Tr(aQW' + (l-a)QW") = 

= a £p/&(Q]M') + (1 - a) ^ P7r Tr(Q«') 

< an t P k + (1 - a)n t Pk 
= n t P k . 

Thus aQ fc + (1 — a)Qfc G Ajjj. an d tne set is convex. 

D. Compactness of the strategy sets -4| IC 

To prove that the strategy sets are compact sets we use the fact that, in finite dimension spaces, a closed and 
bounded set is compact. 

First let us prove that Af lc is a closed set. We define the function g : Af lc — > [0,n t Pk], with 

f(Qk) = £ ftrTr(QW). 

tt&Vk 

We see that g(-) is a continuous function and that its image is a compact and thus closed set. Knowing that 
the continuous inverse image of a closed set is closed, we conclude that Af lc is closed. 

Now we want to prove that the set Af lc is a bounded set. We associate to the tuple of matrices (Q,^) 7T& p K 
the following norm ||Q/t|| = / HQ/t^ll! where ||.||2 is is the spectral norm of a matrix. 



Q^lb = \/ max { A Q(^QMW}r=i- 
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Since for all Q k G -4f IC , Q k is a non-negative, Hermitian matrix we have that: 

ma X {A Qr ,(i)}r =1 <rr(QW)<oo, 

and thus: 



IIQfc^lh = ^/max{A Q u)2(i)}^ =1 = ^/max{A Q w(i) 2 }^ =1 < 
In conclusion the associated norm ||Qfc|| < oo. 



OG. 



Appendix B 

We suppose that there exist two different equilibrium strategy profiles: (Qk,Q-k) G -^fc IC x A-jsP an ^ 
(Q fc ,Q_ fc ) G Al IC x -A^, such that (Qfc,Q_fc) / (Qfc,Q_fc). Then the condition given in the theorem, 
C > is met for the particular choice of (Q' k , Q'_ k ) = (Qfc, Q-fc) and (Qj£, Q" fe ) = (Qfc, Q-jfc). 

By the definition of the Nash Equilibrium, the strategies Qfc, k G /C, are the solutions of the following 
maximization problems: 



max u k (Q k ,Q- k ). 
Q k eAf c 

Thus, Qfc satisfy the following Kuhn-Tucker optimality conditions: 
1) Qfc G Af c , which means that: 



qW = (Q W)H h o 
^^Tr(Q^) < n t P k , 



,Vtt G TV 



'(it) 

2) There exist Afc > 0, and the following Hermitian non-negative matrices of rank 1, , for all ir G Vk, 
such that: 



3) 



A* 



^ P7r Tr(QW)-n 4 Pfc 



Vvr G V K ■ 
V Q ^)Ufc(Qfc,Q_fc) 



,W£V K , 



P7r AfcI - $W 
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Having assumed that (Qfc, Q_fc) is also a Nash Equilibrium, Qfc, with k G /C are the solution of: 

max «fc(Qfc,Q-fc), 

and thus Qfc satisfy the following Kuhn-Tucker optimality conditions: 
4) Q fc G -4fc IC , which means that: 



■kZPk 



5) There exist Afc > 0, k G /C and the following non-negative, Hermitian matrices of rank 1, <& k , for all 
7r G Pk sucn that: 



A/,- 



6) 



tt(8WqW) 

Vvr G P# : 



,Vie^, 



Vqwufc(Qjfc,Q_fc) = p^Afcl-^^' 
Using the third and the sixth optimality conditions, the condition given in ([6]) becomes: 

K 

C = J2 E {p7rAfeTr(Q^ r) ) + P7r AfcTr(Q^) - ^AfcTr(QW) - p^A fe Tr(Q 



TrGPjr fc=l 

Tr(Q<">*<*>; 



< ffx. 

k=l I 

< 0. 



^(Q^*^) + Tc(QW*W) + Tr(Q«$«)} 



^^Tr(QW)-ntP fc 



+ Afc 



From the other four K-T conditions, we obtain that all the terms on the right are negative and thus C < 0. But 
this contradicts the diagonally strict concavity condition and so the Nash Equilibrium is unique. 



Appendix C 

We want to prove that there is no optimality loss when restricting the search for the optimal covariance 
matrices to Q fc G .4f IC such that = WfcPj^Wjf, for all 7r G Vk- Let us consider user k G JC. We have 
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that: 



arg max u k (Qk,Q-k) 



arg max < 



^2 Ptt E 1o S2 



arg max < > p^Elog, 
Q fc G^ SIC ' ^ 



arg max I ^ p^og 2 



MttH 



I + pVH fc Wf W fc HfV H + p ^ VH.Wf Qf } W,Hf V H 



I + pHjfeWfQ^WfcHf + p ^ H,WfQfw,Hf 



arg max < > Elogn 



I + pHfcXj^Hf + p H.WfQ^W.Hf 



(32) 

where we denoted with Xjj^ = Wj^Q^W^. Knowing that the utility function is concave w.r.t. the new 

(n) 

defined matrices Xi ', and the channel matrix has independent entries, we can directly apply the results 
given in E21 to prove that annulling the non-diagonal entries of x! 7 ^ can only increase the values of the 



functions E log 2 



I + pH fe X^Hf + pJ2 H/Wf Qf W,H 



f-7rl 

. In conclusion the optimal matrices Xi are 



diagonal, that we will denote with . The spectral decomposition of the optimal covariance matrices are: 

qP = w*pw wf . 
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